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EFFECTIVE MASS OF ^He ATOM 
PACS 05.30.Jp, 68.35.Rh IN SUPERFLUID AND NORMAL PHASES 


The formula for the temperature dependence of the effective mass of a ^He atom in the su¬ 
perfluid and normal phases is obtained. This expression for the effective mass allows one to 
eliminate infra-red divergences, being applicable at all temperatures, except for a narrow fluc¬ 
tuation region 0.97 < T/Tc < 1. In the high and low temperature limits, as well as in the 
interaetionless limit, the obtained expression reproduces the well known results. The temper¬ 
ature dependence of the heat capacity and the phase transition temperature Tc ^ 2.18 K are 
caleulated, by using the formula obtained for the effective mass. In the framework of the ap¬ 
proach proposed in this work, the small critical index rj is determined in the random phase 
approximation. The obtained value corresponds to the well known result. 
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1. Introduction 

The idea that the transition of liquid "^He into the su¬ 
perfluid state is a manifestation of the Bose-Einstein 
condensation was put forward for the first time by 
F. London [1]. It was the “proximity” of the Bose con¬ 
densation temperature in an ideal gas with helium 
parameters to the transition temperature in real ^He 
that suggested him this idea. Although this interpre¬ 
tation of the phase transition is not free from difH- 
culties [2], it correctly describes, in general, modern 
experiments with cooled gases [3,4]. 

The problems in the theory of liquid ^He, which 
remain unresolved till now, include the calculation of 
corresponding thermodynamic functions in the whole 
temperature interval and the calculation of the tran¬ 
sition temperature into the superfluid state, which 
would agree with the experimental value. At the qual¬ 
itative level, a reduction of the critical temperature 
was substantiated by R. Feynman [5], who introduced 
the concept of effective particle mass. For intuitive 


© I.O. VAKARCHUK, O.I. HRYHORCHAK, 

V.S. PASTUKHOV, R.O. PRYTULA, 2016 

ISSN 2071-0186. Ukr. J. Phys. 2016. Vol. 61, No. 1 
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reasons, he came to a conclusion that, owing to the 
interaction between particles, the effective mass has 
to exceed the atomic one. This conclusion is also valid 
for two-dimensional systems [6]. 

However, there exists a competing mechanism. The 
repulsion at short distances effectively increases the 
system density and, consequently, should increase the 
Bose condensation temperature. This conclusion is 
confirmed by the results of theoretical calculations [7- 
11] and Monte-Carlo simulation carried out for the 
model of weakly non-ideal Bose gas [12, 13[. In or¬ 
der to put the experimental results obtained in the 
“^He-Vycor system in correspondence with the results 
of theoretical calculations, the effect of atomic mass 
renormalization and a shift associated with the re¬ 
pulsive part of the interparticle interaction have to 
be taken into account simultaneously [14]. 

In the literature, the value of effective mass at low 
temperatures was mainly analyzed [15-19]. In works 
[20, 21[, the corresponding temperature dependence 
was obtained within the variational approach. The 
properties of helium in the normal phase were studied 
in works [22,23], where the effective mass of particles 
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was used as a fitting parameter to put the calculated 
structural functions in agreement with experimental 
curves. 

A method of calculation of the effective mass of 
^He atoms in the liquid phase was demonstrated in 
work [24], It allows infra-red divergences typical of 
the phase transition theory to be eliminated. A short¬ 
coming of the proposed approach is the poorly sub¬ 
stantiated extrapolation of the “seed” effective mass 
onto a wide temperature interval (proceeding from 
the corresponding expression obtained for the zero 
temperature) and an incorrect behavior of the ob¬ 
tained effective mass in the critical region. 

Another approach to the calculation of the effec¬ 
tive mass was proposed in work [25]. In its frame¬ 
work, the temperature dependence of the heat capac¬ 
ity was determined. The result obtained turned out 
in much better agreement with experimental data 
than the results of calculations on the basis of a 
“bare” mass. However, the expression for the effec¬ 
tive mass, which was obtained in the framework of 
this approach, does not exclude the mentioned infra¬ 
red divergences, because the ideology of the effective 
mass calculations was not oriented to this purpose. 

This work is aimed at finding such an expression 
for the effective mass, which would eliminate infra¬ 
red divergences and reproduce a correct behavior in a 
vicinity of the critical point (excluding, maybe, a nar¬ 
row fluctuation region). At the same time, it should 
be better substantiated theoretically in a wide tem¬ 
perature region. Another task consisted in obtaining 
the temperature dependence of the heat capacity with 
the use of the new effective mass and in comparing it 
with the previous results. 
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where fh is the effective mass; Sq = /2m, zq, and 

So (q) are the renormalized one-particle spectrum, ac¬ 
tivity, and structure factor, respectively, of the ideal 
Bose gas; Bq = aqSq is the spectrum of elementary 
Bogolyubov excitations; aq = S 2Nvq/{Veq) is 
the Bogolyubov factor, Vq = J dR is the 

Fourier transform of the pairwise interparticle inter¬ 
action potential 4>(i?), 
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is the structure factor of a Bose liquid in the pair 
correlation approximation; and 


Xq — aq th 




( 3 ) 


The distribution of Bose particles with the new spec¬ 
trum looks like 


Zo - l’ 


( 4 ) 


whereas the renormalized one-particle spectrum Sp is 
chosen in the form 


£p — £p /Xp Aq, 


( 5 ) 


2. General Formulas 

While calculating the heat capacity for a many- 
boson system, let us use the expression for its in¬ 
ternal energy in the pair correlation approximation 
[24, 25]. The expression can be obtained by averag¬ 
ing the Hamiltonian with the density matrix found in 
work [26]: 


E = N -h 


1 m 


E 

q/O 


- 1 


dSo{q) 


2 TO ^ 1-h A,S'o((?) 9/3 


where Ap is a correction to the spectrum, which is 
to be determined. The value of Aq depends only on 
the temperature and is actually responsible for the 
activity renormalization. After eliminating infra-red 
divergences, the expression for Ap looks like [24] 

Expression (5) for the renormalized one-particle 
spectrum can also be written in the form 


^ 2m{p) 


( 7 ) 
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where the quantity m{p) is regarded as the total effec¬ 
tive mass of a particle, which depends on the absolute 
value of wave vector p. This form for the spectrum Sp 
was proposed in work [27] in order to exclude infra-red 
divergences. It will be recalled that the effective mass 
TO is formed by many-particle correlations, starting 
from four-particle ones, and, generally speaking, it 
depends on the momentum p. It is clear that we are 
interested in the behavior oi ffi{p) as p ^ 0. As the 
total effective mass, we will understand the quantity 
TO = to( 0). In this connection, let us consider the dif¬ 
ference Ap — Aq as p —> 0 in more details. 

At small p-values, the renormalized spectrum (5) 
can be written in the form [24] 
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work [24]: 
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In expression (12), we isolate the quantity 
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which contains the indicated non-analyticity in whole 
and, simultaneously, is much more convenient for the 
analysis. The next step consists in finding a series ex¬ 
pansion for Aoo in the interval of small p-values and 
confining the series to terms proportional to p^, be¬ 
cause the higher-order terms give no contribution to 
the effective mass owing to equality (11). 

On the right-hand side of equality (13), we change 
from summation to integration: 


In our theory, we use the following expression for the 
temperature dependence of the “seed” effective mass 
TO*, which was obtained in work [25]: 
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where Pq = 2fh/{l3h?) and Pq = Po\/(l — ^o)- In 
the subcritical region, the activity zg = 1. Therefore, 
Pq = 0 here, and the expressions obtained above be¬ 
come a little simpler. 

Let us change the variables: q/p = x and dq = pdx. 
Then 


-k 2 [l/{lieq) - I3eqn{l3eq)[l + n(/3eg)]]|, (10) 

where the notation n{x) = 1/(6"^ — 1) is used. 

It is easy to see that the critical-point divergence on 
the right-hand side of equality (9) originates from the 
integrand at small g-values. This singularity is log¬ 
arithmic, as will be shown later. Such a divergence 
is typical of critical phenomena. Our task consists in 
isolating this singularity and finding a correct expres¬ 
sion for the effective mass. For this purpose, let us 
consider the following equality, which follows from 
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The function Apa;/(1 -I- XpxSoipx)) is finite and tends 
to zero, as a: —>■ oo (at a fixed p). On the other hand, 
the function 


a; Pq/p'^ + (x + 1)^ 
2 p2/p2 + (a: - 1)2 
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Fig. 1. Temperature dependence of the quantity Foo = 
= lim Aoo(p)/£p at p = 0.01. Points correspond to the exact 

expression, and the solid curve to the accepted approximation 



falls down to zero in the interval a; > 1. Moreover, 


where p is the density of the Bose system. On the con¬ 
trary, the form of the series expansion for the struc¬ 
ture factor of the ideal Bose gas depends on the tem¬ 
perature interval. It looks like 
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in the subcritical region (T < T^) and 


So{q) 


rm? 1 



+ o[q) 


(18) 


in the supercritical one. This result follows immedi¬ 
ately from the expression for the structure factor of 
the ideal Bose gas [26] 
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In order to analyze expression (15) analytically, let 
us apply such approximations for the quantities 
and 5o((?), which contain only the expansion terms 
presented above. One can check the adequacy of those 
approximations with the help of the numerical analy¬ 
sis, the results of which are exhibited in Figs. 1 and 2. 

The further analysis of expression (15) will be car¬ 
ried out separately in the sub- and supercritical tem¬ 
perature intervals and at the very critical point, be¬ 
cause the approaches, which should be applied in each 
of those cases, are different. 


P^/p^ + {x + lf 

7 12'' _ 1)2 



dx = 0 


for any Pq- and p-values. Those facts allow us to as¬ 
sert (especially if the matter concerns the critical re¬ 
gion) that only the interval, where the pa;-values in 
the integrand are small, makes a non-disappearing 
contribution to the quantity Aao/p"^ as p ^ 0. 

Now, let us proceed to expanding the quantities 
Xq and So{q) in a vicinity of the zero wave vector, 
q — 0. For Xq, the result is obtained easily: 


3. Calculations in the Subcritical 
Temperature Interval 

In the subcritical temperature interval (T < Tc), ex¬ 
pression (15) reads 
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where 


Xq = f3pvo + o{q), 
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Applying the formula 


where 
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and changing the order of integration, we obtain 
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The denominator of the integrand should be factor¬ 
ized, and the whole integrand should be expanded 
in simple fractions. Then the elementary integration 
over the variable x gives 
1 
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After the corresponding transformations and the in¬ 
tegration over the variable a, we obtain 
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Expanding the obtained expression in p, we find that 
the quantities proportional to p^ originate exclusively 
from the last two terms in the braces. As a result, we 
obtain 
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pI>c ln|xi/x2| 2 

37 r 2 / 3 p(l -I- xr) (xi — X 2 )^ 


When approaching the critical point, one of the 
roots, say X 2 , tends to zero, and we obtain a logarith¬ 
mic divergence for the quantity Aqo in a vicinity of the 
critical point. What is the effective mass in this case? 
Returning to the analysis of the expressions for xi and 
X 2 , we may conclude that there exists a temperature 
Tp, at which xi and X 2 are real-valued quantities. In 
this case, the function arc tanh (p/x 2 ) is no more hnite 
in the temperature interval between Tp and Tc and 
diverges when approaching the critical point. Until 
the quantity X 2 remains complex, the arc tanh func¬ 
tion can be expressed in terms of the trigonometric 
arctangent, which is finite. The temperature Tp can 
be easily found by putting the discriminant of the 
quadratic equation (22) equal to zero. Its numerical 
solution gives Tp « 2.13 K if the critical temperature 
Tc ~ 2.18 K. One can see that this is a very narrow 
interval, which can be interpreted as a fluctuation 
one, i.e. when the fluctuations of the Bose condensate 
becomes comparable with its amount. It can also be 
considered as a region similar to the Ginzburg region, 
where the perturbation calculation method fails. In 
any case, in the framework of our approach, we can¬ 
not draw any proper conclusion about the effective 
mass in this narrow interval. Other methods, e.g., the 
renorm-group approach, are required to analyze this 
region. The numerical analysis testifies that the con¬ 
tribution Aqo to the effective mass is very insignifi¬ 
cant at temperatures below Tp. The analytical form 
of this contribution to the right-hand side of Eq. (11) 
is as follows: 

Po^ ln|xi/x2| . . 

67r2p(l -I- x) (xi — X 2 ) 

4. Calculations at Critical Point 

In order to elucidate the divergence character of the 
quantity Aqo/p^, regarded as a function of p, at 
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the critical point, let us make calculations in this 
case. Let us return to formula (20) and put T = 
in it, which means that no = 0. Then, 




j3p , 
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Again, let us factorize the denominator of the inte¬ 
grand, expand the resulting integrand in simple frac¬ 
tions, and integrate over the variable x. As a result, 
we obtain 
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Aoo gives a non-zero contribution to the one-particle 
spectrum of a Bose liquid at the critical point, we ob¬ 
tain the following equation for the determination of 
7 and p: 
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where xq jp is a root of the equation 
(1 -I- k)xp - 1-7 = 0 , 

and xq = — 7/(1 + k). Changing the variables, 
0 / 13^01 = in ths integral above, we have 
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since 1/ |xo| ~ p (p —t 0). As a result, we obtain 
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where p is the small critical index, and p a character¬ 
istic scale of the wave vector in a vicinity of the criti¬ 
cal point. Taking into account that only the quantity 


The result for the small critical index p was ob¬ 
tained for the hrst time in works [28,29]. The cited au¬ 
thors used a method of expansion in reciprocal powers 
of the order parameter dimensionality. The random- 
phase approximation reproduces only the first term 
of this expansion. Therefore, it is no wonder that the 
result obtained for the small critical index differs from 
the result of Monte-Carlo simulations [30]. 

5. Calculations at Above-Critical 
Temperatures 

At temperatures higher than the critical one, quantity 
(15) acquires the form 


A„o = 


PlzpJi 


qdq 


Hence, we showed that the quantity A^o/p^ di¬ 
verges at the critical point as In [p| (p —> 0). Such a 
singularity is typical of critical phenomena. It can be 
interpreted as a consequence of the expansion of the 
one-particle spectrum of a Bose liquid in a vicinity of 
the critical point: 
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where 7 = 27 / 71 . Differentiating it with respect to p, 
integrating the result again over p, and changing the 
order of integration, we obtain 
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In order to calculate this integral, we symmetrize the 
limits of integration over q, make an analytical con¬ 
tinuation of the integrand into the upper half-plane 
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of the complex g-variable, and close the contour of 
integration by a semicircle of radius R. In the limit 
R -A oo, the integral along the semicircle R equals 
zero, because the power of the integrand’s denomi¬ 
nator is larger by two than the power of the numer¬ 
ator. As a result, our integral is equal to a sum of 
residues at the analytical continuation of the inte¬ 
grand into the upper half-plane times 27rb Only three 
singular points of the integrand fall within this half¬ 
plane: q = p + iRq, q = —p + iPo, and q = iPq. (Note, 
by the way, that the multiplier in the denominator 
with the arctan function does not equal to zero over 
the whole complex g-plane.) As a result of calcula¬ 
tions, we obtain 
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pf^Z[jK TTZ 
47r2/3p 2 p , 


dp\ 
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Without specifying the subsequent rather simple 
transformations, we present the final result for Aqo: 
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We expand the expression obtained in a series in the 
small parameter p and keep only the terms propor¬ 
tional to p^. As a result, we obtain 
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The corresponding contribution to the right-hand side 
of Eq. (11) is as follows: 

p4^o^7 (-97 In'(3) + 8 P 0 + 287 ln( 3 ) - I 67 ) 

27Trp (2Po + 7ln(3))3 ’ 

With the help of the numerical analysis, one can get 
convinced in the smallness of this quantity. Therefore, 
its contribution to the effective mass can also be ne¬ 
glected. 


6. Analytical Expression for Effective Mass 

Taking into account that the quantity Aoo gives an 
insignificant contribution to the effective mass, which 
was demonstrated above, and returning to the calcu¬ 
lation scheme described in work [24], we obtain the 
following expression for the effective mass: 


^ (1 + P(T))’ 


(39) 


where ^ ^ 
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and Vg is the gradient operator. 

Let us expand the operator in a series and 

confine the expansion to first three terms, because 
they give us the required approximation. Making sim¬ 
ple transformations, changing from summation to in¬ 
tegration, and taking the meaning of notations po and 
Pq into account, we obtain the following expression 
for the quantity F{T): 


F{T) 


1 f Xgq'^dq 
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A direct inspection easily verifies that the function 
F{T) equals zero in the limits of both low and 
high temperatures. Therefore, in those limits, fh = 
= m*. Using the results of work [25], we obtain that 
lim m Ki 1.7m and lim m = m. 

T—>0 T —>00 
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7 . Numerical Calculation 
of Effective Mass and Heat Capacity 

Let us illustrate the obtained result in the graphic 
form. For this purpose, we should numerically calcu¬ 
late the ratio rn/m. The corresponding calculation is 
self-consistent, because the expression for m includes 
the quantities So{q), Sq, and Uq, which depend, in 
turn, on m. In practice, this situation implies the ap¬ 
plication of an iteration process, which took 3-4 cy¬ 
cles in our case. 

The calculations were carried out for the equilib¬ 
rium helium density p = 0.02185 A“^, the parti- 



Fig. 3. Temperature dependence of the effective mass of ^He 
atom 



Fig. 4- Temperature dependence of the heat capacity of liq¬ 
uid "^He: (solid curve) theoretical result obtained taking the 
effective mass into account, (circles) experimental data [35—37] 

36 


cle mass m = 4.0026 amu, and the sound velocity 
c = 238.2 m/s in the limit T —0 [32]. The exper¬ 
imentally measured structure factor for liq¬ 

uid "‘He extrapolated to the temperature T —5> 0 [33] 
rather than the Fourier coefficient for the energy of 
pairwise interparticle interaction Vq was used as the 
input information. 

In Fig. 3, the temperature dependence of the effec¬ 
tive mass of ^He atom calculated in the approxima¬ 
tion of pair interparticle correlations is exhibited. On 
its base, using the known formula [34], we also cal¬ 
culated the temperature of the Bose condensation in 
liquid '‘He. The obtained value is Tc w 2.18 K, which 
is very close to the experimental value Tc = 2.168 K. 

While calculating the heat capacity, we used ex¬ 
pression (1) for the internal energy of a many-boson 
system in the pair correlation approximation. We nu¬ 
merically differentiated it with respect to the temper¬ 
ature. Figure 4 demonstrates the temperature depen¬ 
dence of the heat capacity calculated with regard for 
the effective mass. 

8. Conclusions 

An expression for the temperature dependence of the 
effective mass of a ‘‘He atom (in both the normal 
and superfluid phases) is obtained. It allows infra¬ 
red divergences, which are typical of critical phenom¬ 
ena, to be eliminated. The expression for the effective 
mass is applicable at all temperatures, except for a 
narrow fluctuation interval between the temperature 
Tp 2.13 K and the temperature of phase transi¬ 
tion. In the high-temperature limit, as well as when 
the interparticle interaction is “switched-off”, the ef¬ 
fective mass transforms into the “seed” mass of a ‘He 
atom. In the low-temperature limit, we obtain a value 
that coincides with the effective mass of a '‘He impu¬ 
rity atom in liquid ‘He, provided that the “seed” mass 
of a ‘He atom is substituted by the mass of a ‘He one 
[18]. In this context, we note that there is no common 
opinion concerning the effective mass even at the zero 
temperature, to say nothing of a wide temperature in¬ 
terval, because the introduction of this quantity into 
consideration is a phenomenological issue and, to a 
great extent, depends on the approaches applied for 
its calculation [15-17]. 

The behavior of the heat capacity curve theoreti¬ 
cally calculated with regard for the effective mass is 
in much better agreement with the experimental data 
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than if without it, in particular, in the supercritical 
region [31]. In addition, in comparison with the “bare” 
mass, the effective mass obtained in this work gives a 
better agreement with the experimental data for the 
heat capacity in the temperature interval of about 
0.5 K above the phase transition point [25]. 

The application of the effective mass made it pos¬ 
sible to shift the phase transition point from the 
value for the ideal Bose gas to the temperature Tc « 
2.18 K. As was already mentioned, the latter value is 
very close to the experimental one. The “bare” mass 
gives rise to Tc ~ 1.94 K in this case [25]. 

In the framework of the approach proposed in this 
work, we also succeeded in finding the small criti¬ 
cal index ry in the random-phase approximation. The 
obtained value differs rather strongly from the recom¬ 
mended one [30], but simultaneously reproduces the 
well-known result of this approximation [28]. 
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Translated from Ukrainian by 0.1. Voitenko 

I.O. BaKapnyK, O.I. PpuzopnaK, 

B. C. riacmyxoe, P. O. UpumyAa 

E<1)EKTHBHA MAC A ATOM A ^He 
B HAAHJIHHHIH I HOPMAJIbHIH OABAX 

P e 3 K) M e 

3Ha0/i,eHO BHpa3 .gjin TeMnepaTypnoi' aajiejKHOCTi ecjreKTHBHOi 
MacH axoMa ^He b na^gnjiuHniH i HopMajiBnin ((raaax, bkhh p.o- 
BBOjine ycynyTH in^jpauepBOHi poabi^KHOCxi i e aacxocoBHHM 
npn Bcix xeiunepaxypax 3a bhhbxkom BySBKOi (J^jiyKxyapiHHOi 
o6jiacxi 0,97 T[Tc < 1. B rpaHHpi bhcokhx i hhsbkhx xeM- 
nepaxyp, a xaKOxc b rpaHHpi BHKjnoueHHa B3aeMO/i,ii, oxpHMa- 
HHH Bnpaa /i,ae Bi^gOMi pesyjiBxaxn. Ha ocnoBi ec^reKXHBHo'i MacH 
poapaxoBano xig KpnBOi xenjioeMHOCxi, a xaKO>K snangeno xeM- 
nepaxypy (JraaoBoro nepexogy Tc ~ 2,18 K. BuKopncxoByiouH 
aanponoHOBaHHH b po6oxi nigxig, oxpnMaHO anauenHa Mago- 
ro KpnxHUHOro ingexcy Tj b nabgHxceHHi xaoxHunnx 4)a3, axe 
BigxBopioe BjKe BigOMHH peaygBxax gtoro nabgnjKeHHa. 
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